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The transverse spin effects may be helpful to distinguish between scalar 
(J PC = ++ ) or pseudoscalar (J PC = h ) nature of the spin zero (Higgs) 
particle once discovered in future accelerator experiments. The correlations 
can manifest themselves e.g. in the distribution of acollinearity angle of 
in the decay chain H/A — > t + t~; t ± — > i/X . This delicate measurement 
will require however reconstruction of the Higgs boson rest-frame. Then, 
questions of the combined detection-theoretical effects may be critical to 
establish the reliability of the method. An appropriate Monte Carlo pro- 
gram is essential. 

To enable such studies we have extended the standard universal interface, 
of the TAU0LA r-lcpton decay library, to include the complete spin effects for 
r leptons originating from the spin zero particle. The interface is expected 
to work with any Monte Carlo generator providing Higgs boson production, 
and subsequent decay into a pair of r leptons. 

Examples of numerical results and cross checks of the program will be also 
given. In particular, we find that effects of beamstrahlung may be crit- 
ical to the quality of the measurement of the Higgs boson, unless some 
improvements of the method can be found. 
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One of the main goals for future high energy experiments is to measure 
properties of the Standard Model (SM) Higgs sector. Proton-(Anti)Proton 
Colliders, such as Tevatron [1] or LHC [2, 3] are expected to discover the 
Higgs boson, if the SM or one of its AiSSM extensions is true. Otherwise 
the spectrum of possibilities is practically unlimited and discovery can not 
be guaranteed. The comprehensive precise measurements of all Higgs boson 
properties are expected to be left for future experiments on high energy e + e~ 
linear colliders such as JLC [4], NLC [5], or TESLA [6]. 

One of the important measurement, just after establishing that the 
newly discovered particle has indeed spin zero, is to check if it is a scalar or 
pseudoscalar. Depending on the mass of the (to be) discovered Higgs boson, 
and if it is of Standard Model or one of its numerous extensions, different 
observables can give access to this information. Already a long time ago, 
see e.g. [7], it was argued that exploring transverse spin correlations in the 
Higgs boson decay H/A — > can in some cases provide a 

model independent test. The method relies on the properties of the Higgs 
boson Yukawa coupling to the r lepton, which in the general case can be 
written as f(a T + ib T ^^)T (for a discussion of the Higgs boson models, see 
e.g. [8] page 123). The method, at least in principle does not depend on the 
Higgs boson production mechanism at all. 

There are many reasons why this process may turn out not to be in- 
teresting. The cross section may be too small for the luminosity of the 
future collider, the mass of the Higgs boson may be heavy and other Higgs 
boson decay channels better suited for the parity measurement. Finally 
the parity of the Higgs boson may be measurable from the properties of 
its production. However, it is generally accepted that the H/A — > t + t~ 
offers a very interesting signature. Its feasibility needs to be studied espe- 
cially in the context of decisions to be taken on properties of the future LC 
detectors which may be taken soon. The proposed measurement [7] is exper- 
imentally involved. It requires reconstruction of the acollinearity angle (5*) 
between r + r~ decay products in the H/A rest-frame. Note that in case of 
H/A — ► r + r~ the four momentum of the Higgs boson is not directly mea- 
surable as we have the unobservable v T among its decay products; it needs 
to be reconstructed from the constraints of energy momentum conservation 
for the whole event. The distribution in angle (5*) is sensitive to the trans- 
verse t + t~ spin correlations, which are different for the scalar, pseudoscalar 
or the mixed state (we will take into considerations only the extreme cases 
corresponding to choosing either b T or a T equal zero). Precise enough re- 
construction of the H/A rest-frame is important. Many effects, theoretical 
(e.g. QED bremsstrahlung) , or experimental (uncertainty on beam energies, 
not sufficient hermeticity of the detector, angular /energy resolution for all 
particles and jets etc.) may invalidate the method. In the following, we will 
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concentrate on the feasibility of the method, taking into account properties 
of the H/A — > t + t~ decay, and simple assumptions on bremsstrahlung and 
beamstrahlung in reconstructing Higgs boson rest-frame, leaving out from 
the considerations all other limitations and constraints, be it theoretical or 
experimental. 

It is generally expected that the Monte Carlo method is the only way to 
estimate whether the measurement can be realized in practice, and which 
features of the future detection set up may turn out to be crucial. Our 
paper is organized as follows. First, an algorithm for generating decays of 
r ± leptons produced in H/A — > t + t~ including full spin correlations for 
the Higgs boson production mechanism is explained, and some numerical 
examples testing the correctness of the program are given. Later, results 
taking into account inaccuracies in reconstructing the Higgs boson rest- 
frame are shown and conclusions are given. 

Since the Higgs boson spin is zero, the spin correlations of its decay prod- 
ucts do not depend at all on the mechanism of the Higgs boson production. 
Technical difficulties related to the choice of r + and r _ spin quantization 
frames, present in the case of e + e~ — > Z/7 — ► t + t~ [9, 10] (bremsstrahlung 
effects included or not), are not present. The analytical form of the density 
matrix is simple. To calculate the density matrix for the pair of r-leptons it 
is thus enough to: know their four momenta, know that they indeed origi- 
nate from the Higgs boson and, assume the type of the Yukawa interaction. 
Such information is stored in the event data structure called HEPEVT com- 
mon block [11] used by practically all Monte Carlo generators for Higgs 
boson production. 

In Refs. [12,13], the algorithm was developed where all r leptons found 
in the HEPEVT common block can be decayed with the help of the TAUDLA 
library [14-16] and the r decay products are appended to the HEPEVT as 
well. The kinematical information on the momenta of all particles forming 
an event was used to calculate, in some approximation, the longitudinal 
spin state of the r. For our purpose that solution had to be extended, to 
incorporate the full density matrix of the t + t~ pair, in the case when it is 
originating from the Higgs boson decay. The following changes had to be 
introduced to the algorithm explained in Ref. [12]: 

1. The quantization frames for the spin states of r + and t~ need to be 
properly oriented with respect to each other. In our solution they are 
simply connected by the boost along r lepton momenta as defined in 
the Higgs boson rest-frame. At the technical level this is enforced by 
the TRALDR routine [14] defining the relation of the spin quantiza- 
tion frames and the laboratory frame. As an intermediate step this 
routine uses the Higgs boson rest frame. 
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Fig. 1. The ir + ir~ acoplanarity distribution (angle <f>* ) in the Higgs boson rest 
frame. The thick line denotes the case of the scalar Higgs boson and thin line the 
pseudoscalar one. 



2. The density matrix was taken from Ref. [7] and adapted to the quanti- 
zation frames as specified in previous point. Only two cases of purely 
scalar or pseudoscalar Higgs boson were implemented. Any further 
extension is however straightforward. 

3. Generation of r + and r~ decays is then performed following the 
method explained in Ref. [14] and used in KORALB [10] since a long 
time ago. 

4. We have assumed that production generator provides two-body Higgs 
boson decays to r leptons only, in particular, that it does not provide 
any bremsstrahlung corrections. Instead, PHOTOS [17, 18] can be used 
for that purpose, once generation of r ± decays is completed. 

5. More complete inclusion of bremsstrahlung corrections would require 
a substantial re-write and extension of the program to the solution as 
in Ref. [19] or a similar one. 

Once we have explained the main principles of our calculation, let us turn 
to the discussion of numerical results. As an example we will take a Higgs 
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Fig. 2. The n + Tt~ acollinearity distribution (angle 8* ) in the Higgs boson rest frame. 
Full angular range < S* < tt is shown. The thick line denotes the case of the 
scalar Higgs boson and thin line the pseudoscalar one. 



boson of 120 GeV. In the first two plots, which will be constructed for the 
quantities defined in the Higgs boson rest-frame we are totally independent 
of the production mechanism. We will study the predictions for the scalar 
and pseudoscalar cases, essentially to provide the test of our generator. 
Thick lines will denote predictions for the scalar Higgs boson and thin lines 
for the pseudoscalar one. As in Ref. [7], we take the — > i/7r ± decay mode 
only. 

Fig. n presents the distribution in the angle <p* = arccos(n + ■ n~) where 
n 1 * 1 = — i x p — — , i.e. the non-observable acoplanarity angle. The distribu- 

\p n X p T 

2 

tion is indeed, as it should be [7], proportional to ~ iTfg cos</>* respectively 
for scalar and pseudoscalar Higgs. In Fig. [2 we plot the distribution of the 
7r + 7r~ acollinearity angle (5*). The difference between the case of a scalar 
and a pseudoscalar Higgs is clearly visible, especially for acollinearities close 
to 7r (see Fig. 0J). 

Let us now turn to the distributions defined for the semi-realistic case. 
We need thus to take into consideration the combined process of decay 
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Fig. 3. The 7r + 7r~ acollinearity distribution (angle S* ) in the Higgs boson rest frame. 
Parts of the distribution close to the end of the spectrum; S* ~ tt are shown. 
The thick line denotes the case of the scalar Higgs boson and the thin line the 
pseudoscalar one. 



and production of the Higgs boson. As an example 1 , for the production 
mechanism we took the process e + e~ — > ZH; Z — > /j> + fi~(qq);H — > t + t~ 
(only the scalar if can be produced in this process in <SA4), at Center-of- 
Mass-System energy of 350 GeV simulated with PYTHIA 6.1 Monte Carlo 
program [20], effects due to initial state bremsstrahlung were taken into 
account. As in this case production of the pseudoscalar Higgs boson is ex- 
cluded, to quantify the size of the spin effect we will compare the predictions 
when all spin effects are included (thick lines on the following plots), with 
the case when longitudinal spin correlations are included only (thin lines). 
The difference between the two lines visualizes the size of the transverse 
spin effects. 

If we could compare predictions for scalar and pseudoscalar, the differ- 
ence would be roughly a factor of two larger than between the cases of full 



1 We have checked that in case of other production processes and center-of-mass system 
energies, the results, presented later in the paper, remain similar or are slightly less 
sensitive to the transverse spin ( i.e. Higgs boson parity) effects. 
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Fig. 4. The tt + tt~ acollinearity distribution (angle S) in the laboratory frame. Full 
angular range < 5 < tt is shown. The thick line denotes the case when all spin ef- 
fects are included in the decay of the scalar Higgs boson, while only longitudinal spin 
correlations are included for thin line. The two lines are nearly indistinguishable. 



spin, longitudinal spin 2 . Then however, we could not limit our discussion to 
the properties of the Higgs boson decay. Many possibilities due to generally 
distinct, and model dependent, production mechanisms for the scalar and 
pseudoscalar Higgs boson would make the picture more involved and not 
suitable for our discussion. 

As we can see in Fig. |1J the ir + 7r~ acollinearity angle (<5) distribution 
in the laboratory frame looks quite different than in the Higgs boson rest- 
frame, the two cases of different spin treatments are indistinguishable, dis- 
tribution is not peaked at 5 ~ it at all. 

If information on the beam energies and energies of all other observed 
particles (high px initial state bremsstrahlung photons, decay products of 
Z etc.) are taken into considerations the Higgs rest frame can be recon- 
structed. We may define the "reconstructed" Higgs boson momentum as 
the difference of sum of beam energies and momenta of all visible parti- 

2 By numerical accident, the case when only longitudinal spin correlations are included 
is equivalent to the case of non-coherent sum of scalar and pseudoscalar Higgs boson 
contributions of the equal proportions. This holds of course if the same production 
mechanism could be applied for the two cases. 
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Fig. 5. The tt + tt~ acollinearity distribution (angle 5') in the scalar Higgs boson 
reconstructed rest frame. Full angular range < 5* < n is shown. The thick line 
denotes the case when all spin effects are included, while only longitudinal spin 
correlations are taken for thin line. 



cles, that is, for example, decay products of the Z and all radiative photons 
of | cos#| < 0.98. In our study we will mimic in a very crude way beam- 
strahlung effects only, assuming a flat spread over the range of =L 5 GeV for 
the longitudinal component of the Higgs boson momentum with respect to 
the generated one 3 . This assumption means, that the detection effect which 
is practically independent from the Higgs boson production mechanism is 
included only. As we can see (Figs. [5] and EJ) in the distribution of the 
acollinearity angle (5*) defined in reconstructed Higgs boson rest-frame, the 
effects due to transverse spin effects are only barely visible. 

We should stress that, in this very simple example, we have not dis- 
cussed at all other effects potentially degrading the method, such as limited 
statistics, backgrounds, uncertainties in reconstruction of the energies and 
directions for the particles and jets, which may lead to systematic errors 
comparable in size to the parity effect, remaining after beamstrahlung ef- 
fect is taken into account. Alone, this ambiguity in reconstruction of the 
Higgs boson four-momentum degraded the method of measuring the Higgs 



3 The typical spread for the beam energy in linear collider is of the order of few percent 
[4] or even worse. 
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Fig. 6. The tt+tt - acollinearity distribution (angle 8') in the scalar Higgs boson 
reconstructed rest frame. Parts of the distribution close to the end of the spectrum; 
5' ~ 7T are shown. The thick line denotes the case when all spin effects are included, 
while only longitudinal spin correlations are taken for thin line. 



boson parity using the decay chain h — ► t + t~ , r — > tt v in a decisive way. 
We have studied several mechanisms of the Higgs boson productions, in 
all cases depletion of the acollinearity distribution sensitivity to transverse 
spin effect was quite similar. We can conclude that our results are thus 
independent from the production mechanism. 

Recently some work was started on studying detector effects, see [21] 
for details. We can nonetheless conclude that, due to the beamsstrahlung 
effect, there is little hope for the elegant method of reference [7] to check 
Higgs boson parity using its decay to r leptons (whatever the luminosity of 
future Linear Collider), unless, may be, other, unfortunately less sensitive 
to spin than r ± — > -k^v decay modes are used as well. We may hope, also 
that methods similar to the fruitful ones for measurement of r polarization 
at LEP 1, or for the study of CP parity and known since a long time, see 
e.g. [22,23], may become available for our case as well. Definitely, realistic 
studies, combining accelerator, detector and theoretical effects are needed 
to settle the matter. 
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